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SUMMARY 


The minimuiii wave -drag problem with auxiliary conditions is solved 
for axial flow about bodies of revolution consisting of two symmetrical 
ogival sections joined by a circular cylinder. The auxiliary conditions 
are that the total length, the length of the cylinder, the frontal area, 
and the volume are held constant. The results are related to similar 
results known for bodies of revolution without a cylindrical midsection, 
and it is found that the addition of small amounts of center section has 
little effect on the drag. The maximum thickness ratio leading to the 
least total of wave and friction drag is investigated briefly. 


INTRODUCTION 


The formula for the determination of the wave drag of a slender 
body of revolution in a supersonic free stream parallel to the axis of 
the body was first given by von K^man and Moore (reference l) . In a 
later work (reference 2), von Karman reformulated the problem and gave 
the form of the body of prescribed length and max im um cross section 
having a minimum wave drag. The bodies treated in both the above papers 
consisted of ogives at the upstream end of cylinders extending to 
infinity downstream. Somewhat later, Lighthill (reference 3 ) the 

solution to the problem of minimum drag with the auxiliary conditions of 
prescribed length and maximum thickness for a body consisting of two 
symmetrical ogives placed back to back. A paper by Haack (reference 4) 
gives a complete summary of all previous solutions, as well as some new 
solutions, for both symmetrical bodies and bodies of the type discussed 
in references 1 and 2. In reference 5, Busemann has attacked the problem 
of minimum drag of bodies of revolution by exploiting its analogy to the 
problem of wing-trailing-vortex systems of minimum energy. Sears (refer- 
ence 6) discusses the body consisting of two ogives placed back to back 
but, in the case where length and msiximum cross section are prescribed, 
he does not limit the analysis to the case of fore-and-aft symmetry. 
However, the results show that the least drag does occur for symmetrical 
bodies . 
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The present report offers an extension of the results outlined 
above hy taking into consideration bodies consisting of two synmetrical 
ogives joined hy a cylindrical center section,^ (See fig, 1 .) Since the 
stem section is pointed, the question of base drag does not arise. If 
it were desired to consider bodies with finite area at the stem, such 
as boattailed bodies, again without taking into account the base 
pressure, the method used in this report would be applicable, provided 
the meridian section of the body has zero slope at the end >rhere the 
boattall occurs. The analysis then becomes considerably more compli- 
cated than that of the case treated here. 

The introduction of the center section brings a new geometrical 
parameter into the problem, namely, the length of that center section. 
The minimum drag problem can be fomnilated as an isoperimetric problem, 
since the auxiliary conditions are expressible in integral form. It is 
solved under the conditions that frontal area (or maximum thickness), 
volime, length of cylinder, and total length are held constant. This 
rather restrictive set of conditions is then relaxed to include cases in 
which two of the geometric parameters are fixed while the other one is 
free to vary. In this way, three distinct minimum problems connected 
with the type of body considered here can be Investigated systematically. 

Finally, in the appendix, the frictional drag of a body of revolu- 
tion is taken into account in an approximate manner to determine the 
thickness ratio of a body having the least value of combined wave and 
frictional drag. 


LIST OF IMPORTANT SYMBOIS 


B((T,k) i [E(c,k) - k' 2 p(cr,k)l 


Cj)* wave drag coefficient, based on the area 


C])- wave drag coefficient, based on frontal area of body 

f \ 


iThe present work generalizes particu lar cases of bodies with cylindrical 
midsections considered by Max. A. Heaslet and Harvard Lomax in the 
forthcoming series on High-Speed Aerodynamics and Jet Propulsion, 
Princeton University Press. 
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D 

D((T,k) 

E 

E(ff,k) 

F(ff,k) 


k 


k» 


vave drag coefficient, "based on 2/3 power of volume of "body 



(K-E) 
k2 ^ 


i[F(ff,k) - E(c,k)] 


complete elliptic integral of second kind, modulus k 

incomplete elliptic integral of second kind with argument 
a and modulus k 


incomplete elliptic integral of first kind with, argument 
a and modulus k 


modulus of elliptic integrals 



complementary modulus 



K 

21 

2 L 

^o 

r(x) 

R 

So 

So* 

S(x) 

'*'0 


complete elliptic integral of first kind, modulus k 
total length of "body of revolution 

length of cylindrical midsection of "body of revolution 
wiftTiminn radius of "body of revolution 
local radius of "body of revolution 


wave drag divided "by free— stream dynamic pressure 
frontal area of body of revolution (irro^) 



local cross-sectional area of body of revolution [jrr 2 (x)] 


maximum thickness ratio of body (reciprocal of fineness ratio) 
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k 


t 

V 

V* 

Vo 

X 

Z(a,k) 


X,|i 


i 

a 

Po 


local thickness ratio of body 


total volume of body 




velocity of the free stream, 
coordinate along axis of body 

Jacobian Zeta function of argument ct and modxilus k 
Z((T,k) = E((T,k) - I F(d,k) 

Lagrange multipliers (eqiiatlon (4b)) 


X 

I 

argument of elliptic Integrals 



Z^-x2 



free— stream density 


ANALISIS 

Nomenclatxire and Boundary Conditions 


An example of the type of body to be considered in this report is 
shown in figure 1. Also in that figure is shown some of the notation to 
be used. If S(x) denotes the cross-sectional area of the body at any 
point, then 

S(x) = itr2(x) (1) 

where r(x) is the local radliis of the body. It will be stipulated 
that the body is symmetrical about x=0, that it closes at each end, and 
that the ogival sections fair into the cylindrical section with vanish- 
ing slope. Analytically, these conditions become (see fig. l) 

r(±Z) = 0 

r(±L) = ro 

r»(±L) = 0 

where a prime denotes differentiation with respect to x. In terms of 
the sirea fvinctlon S(x), these conditions become (since S*=2nrr’) 
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5(±Z) = 0 
S(±L) = So > 
S' (±L) = 0 


(2a) 


irtiere S© = icco^ is the cross-sectional area of the cylinder, 
cylinder the conditions on S(x) are 


S(x) = So 
S'(x) = 0 I 
S"(x) = 0 


L<x<L 


The cross-sectional area must alvays he positive or zero; 


On the 


(2h) 


S(x) > 0 


(2c) 


Finally, the restriction is made that the maximum cross section occurs 
at the cylindrical portion. Thus, 


S' (x) >0; - I < X < - L 
S' (x) < 0; L< X < I 


(2d) 


The Variational Prohlem 


On the basis of the work of reference 1 or 3# "the wave drag of a 
body such as is illustrated In figure 1 is given by 


R = = 

fpoVo' 



S'(x) dx 



3"(xi) 

X— Xi 


dxi 


In order to arrive at this approximation, it is assumed that the 
body is slender (to «1), and that both S(x) and S' (x) are continu- 
ous and equal to zero at the ends of the body. 


Because of the fore-and-aft symmetry of the body, the above expres- 
sion for wave drag can be modified into one involving integration over 
either the nose or stem section alone, Th\is, for integration over the 
stem (x>0) 




X3-Xi2 


dxi 


(3) 
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The variatioDal prohlem to he solved is of the isoperimetric type, 
since the drag is to he minimized under the auxiliary conditions of 
constant length, frontal area, and volume. The hody shapes determined 
as solutions to this prohlem vill he referred to as optimum bodies. With 
the auxiliary conditions just mentioned, the quantity to he minimized can 
he written 


T — R + V + Sq 


(4a) 


where V is the total volume of the hody, and Xi aud are undetei>- 

mined constants, the so-called Lagrauge multipliers. The volume V can 
he expressed as 


V = 2 


S(x) dx + 2L So = - 2 


X s’ (x) dx 


and the frontal area as 


5q = - y S’ (x) dx 


Equation (4a) can now he written 


T = - 

jt 


r'^ s’(x) 

r' jsiiiiji toi . XX ♦ J 

/l 

7l =‘^1= 


dx 


(h-t) 


where \ and (i have replaced — jtXi and — — respectively. 


In performing a variation of the quantity T, just defined, only 
so-called weak variations will he allowed. This means that the cross- 
section distribution S(x) is deformed slightly, in such a way that the 
derivatives of the deformation function are of the same order of small- 
ness as the deformation function itself. (See reference 7 # P« 7«) The 
variation can he performed in any of a variety of ways, and the resulting 
necessary condition for a minimum (vanishing of the first variation) 
leads to the equation 



2xS”(xi) 

X2-Xi2 


dxi + Xx + p = 0 


(5) 


The function S(x) obtained by solution of this integral equation, and 
two subsequent integrations, is the distribution of cross-sectional area 
which characterizes the required optimum hody. 
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In order to shov that eqmtion ( 5 ) Is the condition for a mininum 
instead of a maximum, the second variation can he examined. It is fovmd 
that the second variation is proportional to the drag of the variation 
of the profile acting alone and, hence, is positive, hy analogy to the 
resiilt found hy Munk in his vork on minimum drag of vings (reference 8 ) . 


Determination of the Cross Section S(x) 


Equation (5) can he written 



2 S* ’ (xi) 


dxi 



( 6 ) 


It is only necessary to solve this equation for x>0 because of the 
symmetry of the body. Make the transformations 


Xi^ = T ; x2 = t 


in equation ( 6 ); it becomes 

s* * (xi) dT ^ _ 

Jjz t-T ~Jr 

Equation ( 6 a) can he written in the form 



( 6 a) 


w(t) 



g(r)dT 
t— T 


which is the familiar airfoil equation. The general solution to the 
airfoil eq\iation is known, being 


g(t) = 




'•(t) y('b-T)(T-g) 

t— T 


dT 


The quantity appearing in the solution as 


g(r)dT 


is of the nature of an arhitrsiry constant, and is to he eval\iated from 
conditions of the problem. In the present case, the constant is 
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P ^ [S(xi)] -^ = 2[S’(Z) - S'(L)] = 0 

by condition (2a) and the condition S'(Z) =0 imposed on equation (3)« 
The solution to eqiiation (6) can he written (in terms of x) 


S' ' (x) 


X 



Z^-X2) (x2^L2) 



The integrations of equation (7) can be performed, yielding 

S..(,) + % f- (x^L^)K-Z^(K-E) . 

2« ^( Z2-x2) (x2-L2) L Z y(Z2-x2)(x2-L2) 


( 8 ) 


where 
K, E 

Z((T,k) 


complete elliptic integrals of first and second kind, 
respectively, modulus k 

Jacobian Zeta function of argument o' and modiiliis k 

Z(a,k) = E(q,k) - ^ F(a,k) 

K 

F(cr,k) E(a,k) incomplete elliptic integrals of first and second kinds, 
respectively, of argument n and modulus k 

k modulus of elliptic integrals 

cr argument of elliptic integrals 


(•■A 


^-X2\ 

z^L2y 


Next, the first derivative, S'(x), cem. be determined by 

S'(x) = / S"(xi) dxi 

It is found that 


( 9 ) 
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Finally, S(x) is given by 


Six) = - 


S'(xi) dX] 


M 

6jt 


2L2F(a,k) - (z2+L2)E(cr,k) + 


I /(i2-x2)(x2-l2) 


+ ^ |kx 2 Z(a,k) + 


Z^E E(cr,k) - F(cr,k) 

This function, S(x), gives the cross-sectional area distribution of an 
optimum body of revolution of the type shown in fig\ire 1, vhen the 
lengths, frontal area, and volume are prescribed. Since the solution 
appeeirs in terms of the undetermined constants X and ^, it is neces- 
sary to find these constants in terms of the prescribed quantities . 

This can be done by determining the frontal area and volume: 


- K Y y (Z^x2)(xa-L2) I (10) 


So = S(L) = ^ [(Z^+L^) E -2L2r ] + (Z^E^-L^K^) 

6it 


(11) 


V = 2L So + 2 



S(x) dx 


= 4 + — 

16 ' 3n 


[ (z2+L^)E -2L^ K ] 


( 12 ) 


Thus, equations (ll) and (12) serve to determine the constants X and 
|i in terms of the prescribed quantities L, 2, V, and So» 


The remaining quantity to be evaluated is the drag. A combination 
of equations (3) and (6) yields 


(13) 

The solution obtained as equation (10) must nov be examined to 
insure that it satisfies the boundary conditions. The conditions of 
equations (2a) have already been imposed in the analysis, as have the 
conditions of equations (2b). The other boundary conditions, (2c) and 
(2d)^ are more complex, however, and require some cstre in application. 
First, notice that if the conditions S(L) S(Z) =0 are met, and 

S’(x) is negative in the interval L £x ^ 2, then certainly S(x) 
cannot become negative in that interval. Thus if the condition (2d) on 


R,i(|v..So 
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the derivative of S(x) is satisfied, then cond.ition (2c) is implicitly 
met. The remaining boundary condition can now be stated as 

S’ (x) = - ^ jJ /( i^-x2)(x2-L2) + ^ KxZ (a,k) < 0 , L < x <Z 


Analysis of this inequality at the end points shows that it iniplies the 
following two conditions on X, and ll: 


XI 

XI 

where 

D 

B 



(l4a) 


(l4b) 

K-E 


k 2 


E-k*®K 



k2 


It is interesting to note the meaning of the two equalities con- 
tained in expressions (l4a) and (l4b) in terms of the body geometry. 
Equation ( 8 ) can be written in the form 



which shows that S’’(x) is infinite at L -unless the equality of 
expression (l4a) holds, and is infinite at i unless the equality of 
( l4b ) holds . Since 

S”(x) = 2n (r'^ + rr»') 


the sing-ularity at L indicates that r’ * is infinite there, while a 
singularity at i indicates that r' is infinite at z. On the other 
hand, if S''(x) is zero at L, then r’ ’ is zero, showing that the 
ogival section fairs into the cylinder -with vanishing second derivative 
as well as vanishing first derivative. Similarly, the vanishing of 
S’*(x) at Z gives a zero value of r’ at the tip, so that the body is 
cusped. Since this only occiurs when the equality 

XI + ^ p = 0 

holds, it is seen that, in general, the optimum bodies have vertical 
tangents at the tips. 

It is convenient to have the formiilas pertinent to the solution for 
optimum bodies in dimensionless terms. Introducing the folio-wing nota- 
tion 


5 



30 = 


So 

Z2^ 


V* 


V 


13 * 


k' 


^ L 


I 
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and 




22 


= Jtt^ 


where t is the local thickness ratio 
put into the following form 


( diameter ) 

V~2l / 


the eqmtions can he 


t3(|) 


B(a,k) - k'^ D (cr,k) 






_y_ 

7t3 


l®KZ(a,k) + E E(ff,k) - k'^KF(cr,k) - K| /fl- 12) (|2_k»2) 


•where | varies from k' to 1, as x varies from L to 2 . 

So*= Jtto2 = ^ k2(B-k'^D) + (E^-k’^K®) 

on it'=^ 

where tQ is the maximum thickness ratio of the body, 

V* = k"* + k2(B-k'^D) 

16 3« 


= 


drag 


1 n V 2?2 ^ 

2 Po''o •' 


= r ( ^ V* + So* 


(16) 

(17) 

(18) 
(19) 


Equations (17) and (18) can he solved for X and |i, resulting in 

(20a) 
(20h) 


X2 = ^ [3(E'^-k'^K^) V* - Jtk^ (B-k'^D) Sq* ] 


where 


^ [3itk^So* - 8 k^(B - k'^D) V*] 

k A 

A = 9(E2-k'^K^) - 8 (B - k’^B)^ 


By use of the res\ilts of equations (20), the inequalities (1^) can he 
expressed in terms of Sq*, V*, and k'j they become 


where 


^ < ^* < £2 
♦1 V* ~ ^2 

<I>3^= 2k2D (B-k'^D) -3 (E^ - k»^K^) 

Sp [3k2D - 4 (B - k'^D] 

$2= 3 (E^ - k'^K^) -2k2 b(B - k’^D) 

i|r - [4 (B - k'^D) -3k^B] 

^ 4 


(21) 
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Figure 2 shows the region defined by the inequalities (21) with 
plotted versus k* . The upper curve represents 



Sq* ^ ^2 
V* '^2 


and the lower, 



h. 


The shaded region between these curves, which will he called the admis- 
sible region, defines the limits within which the parameters k* and 
Sq*/v* must lie in order that the solution for the optimum body satisfy 
the requirement given in expression (2d). That is, for a prescribed 
value of k’, say, the prescribed values of Sq* and V* must he such 
that the ratio Sq*A* ^alls in the shaded region of figure 2. 


Finally, by using equations (20), the formula for the drag coeffi- 
cient Cp* (equation (19)) can he put in terms of the geometric para- 
meters k* , Sq*, V*; 




3(Sq*)^ 

jrk'^A 


2k(E2-k'^K^ -l6jtk^(B-k’^D)^ + Bn^k"* (22) 


DISCUSSION OF SOLUTION 


From the results obtained in the previous section, one can find the 
characteristics of the body of revolution, of the type shown in figure 1, 
having minimum wave drag when the quantities total length (2Z), length 
of cylinder (2L), frontal area (Sq) and volume (V) are fixed. Althovigh 
the semitotal length I no longer appears in the formulas (equa- 
tions (16), (IT), (18), and (22)), having been absorbed into the dimen- 
sionless quantities Cp*, k* , Sq*> V*, it must he remembered that total 
length of the bodies is fixed. It was also found that when all four 
geometric quantities, i, L, Sq^ V, were prescribed, certain limitations 
upon their magnitude must be obse3rved in order to meet bou n d ar y condi- 
tions set forth in expressions (2c) and (2d) . These limitations are 
most simply expressed in terms of the parameters k’ and Sq /V*, where 
k’ is the ratio of the length of the cylindrical section to total length, 
and So*/V* is the ratio of the volume of the cylinder of radius rQ 
and length I to the volume of the body. The permissible range of 
values for k’ and So*A* i^ given in expressions (21), and is shown 
graphically in figure 2. 

Using equation (22), the variation of the drag coefficient Cp* 
with the variables k* and So*/V* can be found. The calculations 
were made for a value of Sq* of jt/lOO, corresponding to a maximum 
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thickness ratio of l/lO. In the accompanying sketch, a three-dimensional 
view of the variation is shown. The two curves in the base plane are 
just those of figure 2, defining 
the admissible region for k' and 
Sq*/Y*. For a given val\ae of the 
length ratio k’ , the drag coeffi- 
cient Cp* varies parabolically 
with Sq*/V*j the minimum occur- 
ring between the extreme admis- 
sible values of So*/v*m With 
increasing k', Cp* increases 
steadily and the rate of increase 
approaches infinity as k' 
approaches unity, A quantitative 
idea of the veiriation is afforded 
by fig\are 3» where the drag coef- 
ficient Cp* is plotted against 
So*/V* for several vailues of the length ratio k' . The curve shown 
for k’=0 agrees with results of reference 4. The increase in drag 
coefficient with k* is seen to be slow for small values of k', indica- 
ting that the greater available volimie restilting from the cylindrical 
center section may warrant acceptance of the slight increase in drag. 

In order to obtain a more detailed view of the variation of drag 
with the geometric parameters, the variational problem that has been 
solved can be reinterpreted. The existence of the limits on the quanti- 
ties Sq*/V* and k’ (expressions (21)) suggests that two of the trio 
k’, Sq*> V* might be fixed, while the third is left free to vary 
within the ascertainable limits determined by the two prescribed values. 
The total length, 21, of the bodies is also fixed, although it appears 
only implicitly. There are three such combinations possible: 

1 . I , k ' , Sq* fixed; V* free 

2. I, V*, k’ fixed; Sq* free 

3 . I, So*, V* fixed; k» free 

These three problems can be stated in physical terms as follows: 

1. Total length, cylinder length, and frontal area fixed; volume free 

2. Total length, volume, and cylinder length fixed; frontal area free 

3 . Total length, frontal area, and volume fixed; cylinder length free 
These three problems will be considered in order. 
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The Three Subcases 


Case 1; Total length, cylinder lengthy and frontal area fixed; 
volume free , - In terms of the dimensionless parameters, this case applies 
to bodies having k’ and Sq* fixed. The admissible range of the volume 
parameter V* can be determined either from expressions (2l) or from 
figure 2. For each value of V* in tbe range so found, an optimum 
body exists, so that a family of optimum bodies is now detennined. The 
variation of drag for members of this family is readily found by use of 
equation (22), or an estimate may be obtained from figure 3« Now some 
one of this family of bodies must give rise to the least drag, and it 
is clear that this member is determined by the condition X equals 
zero.^ The formulas for thickness distribution and drag for the best 
optimum body in this family become (from equations (16) and (19)) 


t^(i) 


"to^ 


^^(cr,k) + E E(a,k) -k’2KF(a,k)] 


Cd^ = 


E^ k’®K^ 




(23) 

(24) 


The volume parameter for the optimum body is given by 


V* = 2^ B-k«^P 

3 E^-k'^K® 


So 


* 


(25) 


It is convenient to have drag coefficients based upon the frontal, 
area and upon the volxmie (to the 2/3 power) of the body, rather than 
i:5)on the ajrea These are, respectively. 



% = r 




i. n V ^ 
2 ^o'^o ^o 


drag 

LPoVo^(V*)^/®Z^ 


lt2 + 2 
_ _ ^O 

IS E^-k*^K^ 

= 2 

It k4(B-k'2D) 


2 


(V*)V3 


If the length L of the center section is allowed to vanish, k’ 
approaches zero, and the last formulas become 


^Dq I = 

^ L=0 

2?he condition that X be zero corresponds to solving the original 
isoperimetric problem of minimum drag with fixed length and frontal 
area, the volume being unspecified. This problem will lead to the 
best body sought for the case 1 under consideration. This result 
could also be obtained from equation (22) by the ordinary method of 
differentiation. Similar remarks apply to the case when |i is taken 
to be zero. 
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% I 

^ L=0 

These results agree with those of references 4 and 6. The drag coeffi- 
cients for the body with a center section can he expressed in terms of 
the above quantities, so that the effect of adding a center section is 
readily seen: 





1 


E^-k»^K 


2 


(26) 


and 


|(V*)-*/3 k4(B-k’2D)2 


(27) 


Figure 4 shows a plot of the quantity Cp against k‘. 

s 

This figure shows that the drag coefficient based on frontal area rises 
slowly when the ratio of length of cylindrical portion to total length 
is small (the drag coefficient having risen 10 percent when the 
cylinder makes up about 10 percent of the body) , but goes up very 
rapidly for bodies on which the cylindrical section makes up more than 
about 50 percent of the body. 


There is another limiting case of some interest for the body with 
k^ and Sq* prescribed, besides the one in which k* vanishes. That is 
the case in which the cylindrical section becomes infinitely long while 
the nose and stern sections have a prescribed length. Thus, both Z 
and L become infinite while (Z-L) remains fixed. The drag resulting 
from such a configuration is given by 


R 


_ drag _ 


8 Sq^ 
Jt(Z-L)2 


which agrees with a result of reference 2 for an ogive of given caliber 
at the end of a semi -infinite cylinder. 


The shape of the best body of a family can be coii5)uted by means of 
equation (23) . Since the thickness t is given by the ratio of diame- 
ter to total length, the ratio of the local radius r to the maximum 
radius Tq is 




g^KZ(g,k) + [EE(g.k) -k»^KF(g,k)] - Kg /(l-l^) (£^-k»^ )' 


E-- - 


k*^ 


1/2 



where 


(28) 
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Plots of this thickness distri- 
bution for several values of k’ 
are shown in figure 5« The 
accompanying sketch shows a plot 
of eq^uation ( 25 ) in relation to 
the extreme admissible values of 
Sq*/V* versus k’ . The upper 
c\arve corresponds to the eqioality 
in expression (l^b), and repre- 
sents the relation for a body 
with cusped tips, while the low- 
est curve corresponds to the 
equality in ( 1 ^+a), and gives 
zero curvatiire 'v^ere the end 
sections join to the cylinder. 
Therefore, the best bodies of 
the present families never have 
cusped tips, nor zero c\irvature 
where the ogival sections and 
center section join. 

Case 2; Total length, 
volume, and cylinder length 
fixed; frontal area free .- In 
this case, which can be analyzed 
in the same way as the one just 
preceding, the volume and length- 
ratio parameters V* and k’ 
serve to determine a range of 
permissible values of the frontal area parameter Sq*« This range is 
again obtainable either from expressions (21), or from figure 2. Fig- 
ure ( 3 ), showing variation of drag with k’ and Sq*/V*, can once more 
be consulted for a general view of the behavior of the drag as the 
parameters change. 



The optimum body in a given family is now determined by the vanish- 
ing of |i (see footnote 2 ), and the eqTiations for thickness distribution 
and drag become 

t 2 (|) = hi |k2^(a,k) -k»2D(a,k)] -I /(l-|2)(|^k’2)| ( 29 ) 

k2(B-k’^D) I J 


Up - 


3^1 tpav* 
k^(B-k'^D) 


The relation between Sq* and V* for the best member of a family is 


(30) 
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Again the drag coefficient can he based either v5)on (volume)^/^ or 


upon frontal area, the former being the more useful in this case since 
volume is fixed. The drag coefficient based on (volume)^/® is 


In both cases, the function of k* in sqiiare brackets reduces to unity 
as k' vanishes, and the resulting expressions for bodies without 
center sections agree with res\ilts of references and 6. Figvire 6 
shows a plot of the q\aantity 


as a function of the length ratio k’ . The behavior is qualitatively 
the same as fo\md in case 1 for the drag coefficient based on frontal 
area, but the increase of drag coefficient with k’ is somewhat slower 
in the present case. 

In the present instance, where V* and k' are given, the thickness 
distribution function is, from equation (29), 



(32) 


and that based on the frontal area is 



1 1 
(B-k*^D)^J 


( 33 ) 



k4 
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that 


Plots of this thickness distribu- 
tion are shown in figure 7> where 
shapes for several values of k’ 
are given. The sketch shows in 
this case that again the best body 
of a family is one having neither 
cusped nose and stem sections nor 
zero curvature at the junctions of 
end and center sections; the best 
body at k' = 0 has zero curvature 
at the junct\rre of the ogival sec- 
tions . 

Case 3 : Total length, frontal 

area, and volume fixed; cylinder 
length free .- In this case the 
dimensionless parameters Sq* and 
V* are fixed, while the length 
ratio k’ is free to vary. It 
is the most difficult case to 
analyze because the parameter k’ 
appears implicitly in the func- 
tions involved. By differentia- 
ting equation (22) with respect 
to k', holding Sq* and V* fixed, 
it is found that the derivative 
8Cp*/Sk* is never negative, and 
vanishes for a value of k* such 


Sq* 


where $1 and are defined after expressions ( 21 ). In terms of the 
admissible range of the length ratio k* for fixed values of the frontal - 
area and volume paxameters Sq* and V*, this equation means that the 
member of the family for which k' has the least admissible value 
(fig. 2) is the one with least wave drag in that family. It can be seen 
from figure 2 , however, that if the ratio of the given parameters Sq*/v* 
is greater than 8/3« = 0 . 849 , the least admissible value of the length 
ratio is always k’ equal to zero. 


In the cases where the ratio Sq*/V* is less than 8/3«> the drag 
coefficient of the best body of a family is given by 


V = 


(35) 
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vhere the val\ie of k» to be used is the one which makes 

Y* 


( 36 ) 


This vaJ.ue of k' is best fovin'd from figure 2. It is not so useful in 
this case to refer the drag coefficients to those for a body without 
center section, for the value of the length ratio k' is no longer arbi- 
trary and hence cannot be made to vanish at will. 

For the cases in uhich the ratio of fronteil-area parameter to 
voliJine parameter is greater than 8/3jt, the best body of a family is the 
one for which the length of the cylindrical section is zero, as has been 
noted. The drag coefficient Cjj* is then given by the formula 


{2k Y*^ - I6it V* So* + 3^ So*2) (37) 

Since the length of the cylindrical center section is zero and. no longer 
enters as a parameter^ it is convenient to reintroduce the semi total 
length I into the formulas • The wave drag, divided hy free -stream 
dynamic pressure, is 

(2^ - l6n I VSo + (38) 

and this agrees with the result of reference If. The length I must be 
between the limits 


3^ JL <i< It X 

3« So « So 


(39) 


Since for a given cross-sectional area, the body of least wave drag will 
be the one with the longest admissible length, it is clear that for the 
present case, where So and V are prescribed, and no center section 
exists, the best body is the one for which 


= it X 

So 


(ko) 


as stated in reference k. Bodies, the length of which is greater than 
the value of equation (4o), do not fall within the admissible region and 
hence violate condition (2c). 


Shapes of the bodies for the present case are shown in figure 8. 

The body for which Sq*/y* equals k/j: and k' = 0 is shown at the top of 
the fig\n:e. This body is the best one of all those without a midsection, 
having the relation between Z, Sq, and V given in eqtiation (to). The 
other body shapes shown are each the best body for the prescribed value 
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of the parameter ratio Sq*/y*, the value of k' being chosen from the 
lower curve of figure 2. The equations for the shapes are 


for 0.5 


Sq* < 8 
Y * 3jt 



3(E2_k»%2^ 


8k^(B-k*^D) I (E^-k’^K^) ^ - k^(B-k'^D)J^Z.^% 


_ 8 k^(B-k' D) Vi. 

« So* 


2 .1/2 


(4la) 


where (r/ro)j^ and (r/rQ) 
respectively, ^ 

for 



8 ( 3 


V* 

So* 


-1 


are defined in equations ( 28 ) and (3^) 


8 < 8 q * <4 
3« V* ^ 






(to) 


In order to obtain the best body of tbe family described by equa- 
tion (4la), the value of k* to be used can be found from equation ( 36 ), 
or from the lower curve of figure 2. For the family described by equa- 
tion (^Ib), the best member is the one for which 

So* , If 

V* 3 t 


SUMMARY OF RESULTS 


For convenience, the important drag formulas of the preceding 
analysis have been gathered together in the present section. The equa- 
tions are numbered just as they appear in the text. The formulas are 
given in terms of the dimensionless parameters k^, Sq*^ V*, which 
are related to the total length 2l, the cylinder length 2L, the frontal 
area Sq, and the volume V by means of the equations 


k» 

So* 

V* 


Ji. 

I 

12 

X. 
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Case 1: Total Length, Cylinder Length, and 

Frontal Area Fixed; Volume Free 


The cylinder length and frontal area determine the dimensionless 
parameters and Sq*^ vhile the volume is represented by the para- 
meter V*. The prescribed values of k’ and Sq* determine an admissible 
range of values of V* (see fig, 2), and the optimum body of this family 
has the folloving characteristics: 


V* = 2^ So* 

3 E2-k*^K^ 


Drag coefficient based on frontal area, 

Ct) 




(See fig. 


(25) 


(26) 


Drag coefficient based on (volume)^/® 

2 (V*)*/3 k'*(B-k'^D)^ 

JT 


(27) 


Case 2: Total Length, Volume, and Cylinder 

Length Fixed; Frontal Area Free 


The given volume and cylinder length determine the dimensionless 
parameters V* and k^, and these values lead to a range of admissible 
values for Sq*. (See fig. 2.) The best body of the family so deter- 
mined has the following characteristics: 


So* - A B-k»^P 


Drag coefficient based on (volinne) 


2/3 

) 


(31) 


I 


= (See fig. 


6 ) 


(32) 


Drag coefficient based on frontal area. 


9 2 

8 ° 


(B-k’^D)‘ 


( 33 ) 
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Case 3: Total Length, Frontal Area, and Volume 

Fixed; Cylinder Length Free 


The frontal area and volume determine Sq* and V*, and these in 
txjrn determine a range of admissible values of k’,(See fig. 2.) If 
the quotient Sq*/v* is less than 8 / 3 « =0.8i<-9^ the best body is one 
with a value of k' such that 


<S>i Sq* 
Vi V* 


where $1 and \|^i are defined in expressions (21). This value of k’ 
is most easily found from figure 2. The drag coefficient for the opti- 
mem body is then 


Cd* 


nk'^'A 


2l4.(E2-k'^K^)(V*)^-l6jtk^(B-k»^D)V* Sq* + 


( 35 ) 


In case the value of the ratio Sq*/V* is such that 

_8_ < So* < 4 

3:t “V^ It 

the best body is one with no center section (k* = O), and the drag 
coefficient is 


or. 


in terms 
R 


2U(V*)^-l6it V* Sq* + 3«^(Sq*)^ 

of the remaining three parameters I, Sq^ V, 

= ^ ( 21 ^ v2-l6tr I V So + 3«2 z^Sq^) 

2 Po^o 



( 37 ) 

(38) 
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APPENDIX 

A MINIMUM DRAG PROBLEM WHICH INCLUDES AN APPROXIMATION 
TO TEE EFFECT OF SKIN FRICTION 


It is possible to include the effects of skin friction in the 
analysis of minimum drag of bodies of revolution, provided the s\irface 
area of the bodies is known. If it be assumed that this surface area 
is expressible in the form 

Area = 2L(2jtro) + 2(1 -L) (27ro) (Al) 


where the second term represents the area of the ogival end sections, and 
7 is an unspecified constant, then the drag dvie to skin friction is 

=‘*r[l+ (f -1 ) (a 2) 

where Cp^ is the friction-drag coefficient, to is the maxim\mi thick- 
ness ratio, and k' is again the ratio of length of cylindrical section 
to total length. 


Consider now a body of revolution with prescribed length ratio k* 
and frontal-area parameter Sq** The wave-drag coefficient based on 
frontal area for the best such body is, from equation (26) of the text. 



to2 


The total-drag coefficient based on frontal area, is then 


p = "^0^ . V 

^Drj. E^-k'^K^ « "tc 


1 + ( J -1 Ik 




(A3) 


Now the total drag can be minimized with respect to maximum thickness 
ratio to> there results for the optimum, to*, 

1/3 


to' = ^ { 


16 7 Cp^ (E^-k’V) 


1 + f ^ -1 ] k 




(Ak) 


It remains to assign values to the constants 7 and C-^ . 


Inspection of figure 5 shows that the bodies under consideration do 
not differ greatly in shape from prolate ellipsoids of revolution. Thus, 
for the present purposes, it will be sufficient to use the approximation 

-I 


2 k 
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corresponding to the surface area of a prolate ellipsoid the minor axis 
of which is small compared to its major axis. Equation (a 4 ) for the 
optimum thickness ratio now becomes 

^o’ = ^ (E 2 -k>^K^) [1 + (|«- 1 ) k»]| (k 5 ) 

If the cylindrical section is allowed to vanish, k’ vanishes, and the 
last eq\iation becomes 

*°'Jo ’ (W) 

The optimuiDL thickness ratio for bodies with a center section can there- 
fore be expressed as 


to* = c(k*) to* ] 
Jj— 0 


(^ 7 ) 


where 


c (k*) 




A graph of c(k’) versus k* is shown in the sketch. 



The remaining constant is the 
skin-friction drag coefficient 

For the purpose of illustra- 
tion, an average value of 0.0025, 
corresponding to a turbulent 
boundary layer at a Mach number of 
about 1.7 and a Reynolds number of 
13 million, was taken from the 
data of reference 9- The optimiom 
thickness ratio for a body of 
revolution with no center section, 
and with prescribed length, con- 
sidering both wave and friction 
drag, is then found to be (from 
equation (a6) ) 


1^0 


n* 1 

= 1^ 


By means of the sketch of the variation of c(k'), the resTilts can be 
extended to bodies of revolution with a center section. Consider a body 
the center section of which makes up 10 percent of its total length. 

From the sketch and the above value of to* | , it is seen that the 

L=0 
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optimum value of thickness ratio remains about l/l^, again showing that 
the effect of the added cylindrical portion upon the drag is small for 
small values of the length ratio k' . The optimum thickness ratio 
decreases to I/19 for a body the center part of which is 50 percent of 
the total length. 
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Figure / — Type of body considered and notation used. 
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Figure 2- Region of admissible values for the parameters k' and So/V* 



Figure 3 -Variation of drag coefficient with S^/V* for several values of k'. So‘ ir/tOO. 
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Figure 4.^ Variation of drag coefficient based on frontal 
area with length of cylindrical sect ion, given length - 
ratio and frontal-- area parameter. 
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Figure 5- Variation of body shape with k; given length- ratio and frontoi-area pa ram eter. 
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Figure 6.— Variation of drag coefficient based on volume 
with length of cylindrical section, given length-^ 
ratio and volume para m et e r. 



34 


NACA TN 2535 




V? 


NACA - Lanelev Field Va 


Figure 8— Variation of body shape with So/V* given volume and frontal-area parameters. 


